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ABSTRACT 


High-resolution  seismic  velocity  and  attenuation  models  of  the  Tibetan  Plateau  and  adjacent  regions  of  western 
China  are  critical  to  monitoring  seismology  efforts  in  Eurasia.  Our  waveform  data  are  collected  from  new  deployed 
seismic  networks  including  NETS  (North-Eastern  Tibet  Seismic  experiment),  INDEPTH-IV- ASCENT, 
INDEPTH-IV-UK,  and  previous  networks,  comprising  an  unprecedented  coverage  of  the  eastern  half  of  the  Tibetan 
plateau.  Three  seismic  methods  are  used  in  the  velocity  structure  in  this  region:  Pn  velocity  tomography,  body-wave 
finite-frequency  tomography,  and  Rayleigh  wave  phase  velocity.  We  observe  consistent  anomalies  in  the  crust  and 
upper  mantle  using  all  of  these  methods.  These  anomalies  include  a  low  velocity  in  northern  Qiangtang  and  the 
Songpan-Ganzi  fold  belt  and  the  high  velocity  in  the  Qaidam  Basin.  QLg  and  QPg  models  have  been  determined  using 
a  Reverse  Two-station/event  Method,  which  shows  a  high  seismic  attenuation  zone  along  the  Kunlun  belt.  We  have 
also  observed  significant  azimuthal  anisotropy  of  seismic  Q.  The  isotropic  seismic  Q,  suggested  as  an  estimate  of  Qs 
for  the  crust,  is  low  in  northern  Qiangtang  and  the  Songpan-Ganzi.  The  high-Q  directions  parallel  major  fault  planes 
in  eastern  Tibet,  and  correlate  with  fast-direction  of  shear  wave  splitting  and  Rayleigh  wave  azimuthal  anisotropy 
fast  directions,  which  suggests  a  coherency  of  deformation  between  the  crust  and  upper  mantle  in  this  region. 
Similarly,  a  high  velocity  and  high  Q  block  in  southeastern  Tibet  around  eastern  Bangong-Nujiang  Suture  and 
Eastern  Himalaya  Syntaxis  correlates  well  with  the  high  velocity  in  this  region. 
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OBJECTIVES 


The  objectives  of  this  project  are  to  establish  high-resolution  three-dimensional  seismic  velocity  models 
and  regional  wave  attenuation  models  for  the  crust  and  upper  mantle  in  Tibetan  Plateau  and  adjacent 
regions  of  western  China.  The  Tibetan  Plateau  and  adjacent  areas  have  been  focused  by  many  studies  of 
seismic  velocity  structure  (e.g.,  Kind  et  ah,  2002;  Tilmann  et  al.,  2003;  Li  et  ah,  2008)  and  seismic 
attenuation  structure  (e.g.,  McNamara  et  al.,  1996;  Phillips  et  al.,  2000;  Xie,  2002).  Reliable  high  resolution 
images  of  the  crustal  and  mantle  of  this  region  is  very  important  but  limited  by  insufficient  data  sets  and 
mutually  independent  techniques.  Our  research  involves  new  waveform  database,  which  have  increased  the 
coverage  of  the  northern  Tibet  and  Qaidam  Basin  (Figure  1).  All  stations  are  equipped  with  broadband 
CMG-3T,  CMG-ESP,  and  STS-2  sensors.  Based  on  the  results  of  last  year’s  research  review  (Sandvol  et 
al.,  2010),  we  have  applied  this  new  data  set  and  various  techniques  to  estimate  the  seismic  velocity  and  Q 
variations  within  this  area. 


Figure  1.  A  map  showing  the  broadband  stations  from  all  permanent  and  temporary  networks  whose 
data  are  used  in  our  project.  Inverse  triangles  represent  the  experiment  ASCENT  and  stars 
represent  NETS.  As  other  temporary  networks  deployed  across  the  Tibet,  the  stations  of 
Namche  Barwa,  MIT-China,  INDEPTH-IV-UK  are  represented  by  triangles,  squares,  and 
circles,  respectively.  The  black  diamonds  represent  the  permanent  network  CDSN. 
Different  colors  are  used  to  illustrate  the  details  of  seismometers  and  dataloggers.  The  red, 
green,  blue,  purple,  and  yellow  symbols  represent  Streckeisen  STS-2,  Guralp-CMG3T, 
Guralp-CMG3ESP,  Guralp-CMG40T,  and  Guralp-CMG3TD  (only  in  NETS  experiment) 
seismometers  associated  with  Reftek-130  dataloggers;  the  orange  symbols  (only  in  NETS 
experiment)  represent  Guralp-CMG3TD  seismometers  associated  with  Guralp  dataloggers; 
the  small  purple  triangles  (only  in  Namche  Barwa  experiment)  represent  stations  where 
Guralp-CMG40T  seismometers  replaced  other  seismometer  types  in  deployment. 

RESEARCH  ACCOMPLISHED 


Using  both  travel-time  body  and  surface  wave  tomography,  we  have  created  a  Pn  velocity  model,  3-D  P 
and  S  wave  velocity  models,  and  Rayleigh  waves  phase  velocity  models  for  the  crust  and  upper  mantle  of 
this  region.  We  have  also  generated  2-D  QLg  and  QPg  models  while  site  and  source  responses  are  isolated 
and  resolved  as  well.  A  significant  azimuthal  anisotropy  of  Q  is  dramatically  observed  and  solved  in  2-D 
tomography. 
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Three  Velocity  Models 

Pn  waves  mainly  travel  within  a  layer  closely  beneath  the  Moho  discontinuity,  and  thus  its  velocity  model 
mainly  represents  the  velocity  structure  of  the  uppermost  mantle.  We  have  applied  the  method  of  Hearn  et 
al.  (1991)  to  generate  the  Pn  velocity  tomographic  model  (Figure  2).  The  Pn  velocity  varies  from  7.8  km/s 
to  8.3  km/s  with  an  average  of  8.1  km/s.  In  the  Qiangtang  and  Songpan-Ganzi  (Bayan  Har)  terrane  of 
northeastern  Tibet,  low  Pn  velocities  are  generally  observed.  A  high  Pn  velocity  anomaly  is  located  at  the 
eastern  end  of  the  Bangong-Nujiang  Suture  around  the  Eastern  Himalayan  Syntaxis.  The  Qaidam  Basin  and 
Gonghe  Basin  have  high  Pn  velocities.  An  abrupt  Moho  depth  change  is  suggested  by  the  observed 
significant  difference  of  station  delays  along  Kunlun,  the  northern  margin  of  the  plateau. 


Figure  2.  Pn  wave  velocity  structure  in  the  northeastern  Tibetan  Plateau.  High  Pn  velocity  anomalies 
are  located  to  the  north  of  the  Kunlun  range  and  around  the  Eastern  Himalayan  Syntaxis. 
Low  Pn  velocity  anomalies  are  mainly  observed  in  central  Qiangtang  and  northern 
Songpan-Ganzi  terrane. 

We  have  applied  a  finite-frequency  tomography  method  (e.g.,  Hung  et  al.,  2004;  Liang  et  al.,  201 1)  to 
generate  three-dimensional  P  and  S  body  wave  velocity  models  for  the  lithosphere  of  the  eastern  Tibet 
(Figure  3).  A  contrast  between  a  widespread  slow  velocity  northern  Tibet  and  heterogeneous  southern  Tibet 
is  significant  in  the  upper  mantle  velocity  structure.  The  homogeneous  low  P-  and  S-  velocity  anomalies  in 
the  crust  and  upper  mantle  of  the  northern  Qiangtang  and  Songpan-Ganzi  terranes  probably  penetrate  into 
the  upper  mantle.  High  P-  and  S-wave  velocity  anomalies  are  generally  distributed  within  the  mantle 
beneath  the  Qaidam  Basin  and  southern  Tibet  to  a  depth  of  more  than  250  km.  A  few  N-S  oriented  low 
velocity  bands  are  observed  beneath  southern  Tibet  and  appear  to  extend  to  at  least  150  km  depth. 

A  three-dimensional  model  of  fundamental  mode  Rayleigh  waves  phase  and  shear-wave  velocities  with 
azimuthal  fast  directions  is  also  presented  for  the  lithosphere  of  the  eastern  Tibet  (Figure  4).  Like  Pn  and 
body  wave  results,  low  velocity  anomalies  occur  across  and  within  major  strike-slip  fault  zones  in  the 
Qiangtang  and  Songpan-Ganzi  portion  of  the  northern  Tibet.  A  high  velocity  anomaly  is  observed  along  the 
Bangong-Nujiang  Suture  at  65  km  depth,  probably  is  a  portion  of  the  underthrusting  Indian  plate. 

Azimuthal  fast  directions  are  consistent  at  all  depths  up  to  approximately  200  km,  which  suggests  a  vertical 
coherent  deformation  in  the  Tibetan  crust  and  lithospheric  mantle.  The  highest  magnitude  of  azimuthal 
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anisotropy  is  observed  in  northern  Tibet  with  a  direction  nearly  paralleling  major  fault  trend  and  shear 
wave  splitting  fast  directions. 


Vp  model 


-2.50-1.25  0.00  1.25  2.50 


Vp  (%) 


Vs  model 


Vs  (%) 


Figure  3.  Three-dimensional  P-  and  S-wave  velocity  models  for  the  upper  mantle  of  the  eastern  Tibet.  A 
significant  N-S  contrast  for  both  P-wave  and  S-wave  velocity  is  significant. 


All  of  the  methods  individual  results  consistently  show  a  low  velocity  zone  in  the  northern  Qiangtang  and 
Songpan-Ganzi  terranes,  and  relatively  high  velocity  in  southern  Tibet  and  Qaidam  Basin.  The  strong 
correlation  between  the  independent  methods  suggests  that,  at  least  to  first  order,  these  models  are 
relatively  robust.  Further  studies  are  expected  to  see  if  these  models  will  help  to  improve  sparse  network 
earthquake  locations. 


Seismic  Attenuation 

QLg  is  typically  assumed  as  an  approximation  for  shear  wave  Q  of  the  crust,  and  is  widely  used  in  nuclear 
discriminants.  We  used  a  Reverse  Two-station/event  Method  (RTM)  (Chun  et  al.,  1987)  to  measure 
accurate  path-based  1/QLg  values.  This  method  is  an  improvement  upon  the  Two-Station  Method  (TSM) 
and  includes  two  specific  cases:  Reverse  Two  Station  (RTS)  path  and  Reverse  Two  Event  (RTE)  path.  In 
theory,  it  eliminates  the  contribution  on  spectral  amplitude  from  the  source  excitation  function,  radiation 
pattern,  instrument  response,  and  site  response,  and  has  been  presented  as  the  most  accurate  method  in  Q 
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measurement  (Ford  et  al.,  2008).  Because  of  its  rigorous  geometrical  requirement  and  the  specific 
propagation  characteristics  of  Lg,  the  application  of  this  method  still  requires  extensive  seismic  network 
deployed  contemporaneously  throughout  Tibet  (Fan  and  Lay,  2003).  The  deployments  of  ASCENT  and 
NETS  have  ideally  increased  the  accuracy  and  resolution  of  1/QLg  measurements,  and  furthermore,  permit 
an  application  of  both  the  RTS  and  RTE  in  Tibet. 


Figure  4.  The  3D  model  of  Rayleigh  waves  phase  velocities  with  azimuthal  anisotropy.  Notice  the 
color  scale  of  each  section  represents  different  velocity  range.  Low  velocity  anomalies  and 
high  azimuthal  anisotropy  magnitude  mainly  locate  within  northern  Tibet. 
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8 


normalized  log(l/s) 


Figure  5.  The  nearly  linear  relationship  between  logarithm  of  the  reciprocal  of  normalized  RTM  Lg 
source  term  and  magnitude  ML.  The  solid  line  represents  the  least-square  linear  fit  line 
and  the  two  dashed  lines  parallel  to  it  represent  its  95%  confidence  interval. 

We  also  use  the  RTM  to  solve  the  site  response  and  source  response.  A  significant  site  response  difference 
is  observed  between  the  Tibetan  Plateau  and  major  basins,  probably  dependent  on  topography  or 
sedimentary  thickness.  This  is  very  important  to  studies  on  moment  and  ground-motion.  Also,  a  nearly 
linear  relationship  is  observed  between  the  logarithm  of  the  reciprocal  of  normalized  source  term,  solved 
from  RTM  and  Lg  spectra,  and  magnitude  ML  between  3  and  7  (Figure  5).  The  seismic  moment  M0  can  be 
estimated  from  ML  (Hanks  and  Kanamori,  1979).  Because  the  RTM  eliminates  the  contribution  of  path, 
site,  and  radiation  pattern,  it  is  applicable  in  accurately  detecting  the  source  magnitude. 
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Figure  6.  Imaged  variations  of  isotropic  Q,  anisotropy  magnitude,  and  high-Q  direction  in  the  crust 
of  eastern  Tibetan  Plateau  at  1  Hz,  represented  by  the  color  scale  and  the  length  and 
orientation  of  black  vectors,  respectively.  Major  faults  and  sutures  are  also  shown,  while 
strike-slip  faults  are  green,  other  types  of  faults  are  gray,  and  sutures  are  yellow. 

We  observe  that  a  strong  sinusoidal  (or  cosinusoidal  with  a  phase  shift)  relationship  between  the  azimuth 
and  RTM  1/QLg  values  for  each  inter-station  or  inter-event  paths.  Such  azimuthal  variation  in  1/QLg  has 
been  observed  in  the  western  United  States  (Sutton  et  ah,  1967;  Phillips  et  ah,  2009),  Tien  Shan  (Martynov 
et  ah,  1999),  and  southern  Tibet  (Reese  et  ah,  1999).  Similar  correlations  were  derived  and  widely  used  for 
body  and  surface  waves  (e.g.,  Backus,  1965).  We  have  generated  a  2-dimensional  tomographic  model  of 
1/QLg  with  an  azimuthal  anisotropy  by  solving  the  inverse  problem  similar  with  Pn  and  surface  wave 
azimuthal  anisotropy  tomography  (e.g.,  Hearn,  1996).  Approximately  48,000  four-spectral  amplitude 
groups  are  collected  as  input  data  for  the  tomography.  Resolution  tests  suggest  that  an  anomaly  with  the 
size  of  2°x2°  can  be  retrieved  for  the  eastern  Tibet.  Figure  6  shows  the  model  at  a  frequency  of  1  Hz. 

The  1/QLg  azimuthal  anisotropy  is  weakly  frequency-dependent  within  eastern  Tibet.  There  is  a  strong 
azimuthal  anisotropy  of  l/QLg  within  the  northern  Tibet.  The  maximum  anisotropic  magnitude  is  observed 
to  be  higher  than  0.01  in  the  northwestern  Songpan-Ganzi  terrane,  Kunlun  range,  and  northern  Qiangtang 
terrane,  while  the  anisotropic  magnitude  is  very  low  in  the  central  Qaidam  Basin,  showing  a  significant 
discontinuity  along  the  northern  margin  of  Tibetan  Plateau.  The  high-Q  direction  is  nearly  uniformly  E-W 
in  the  northern  Tibet,  very  similar  to  the  pattern  of  surface  wave  azimuthal  anisotropy.  Furthermore,  the 
high-Q  and  fast  directions  parallel  major  strike-slip  fault  planes,  such  as  the  North  Kunlun  Fault  and  South 
Kunlun  Fault.  The  anisotropic  magnitude  increases  around  the  Eastern  Himalayan  Syntaxis,  where  the 
high-Q  direction  likely  has  a  complicated  rotation  from  E-W  in  the  west  to  NW-SE  to  the  east.  By  isolating 
the  anisotropy  of  l/QLg,  the  isotropic  l/QLg  is  strongly  frequency-dependent  in  the  eastern  Tibet.  The 
Qaidam  Basin  has  a  stable  low  isotropic  l/QLg  (<1/500  at  1  Hz)  which  indicates  a  cold,  stable,  and 
relatively  homogenous  crust.  In  northern  Tibet,  such  as  the  Songpan-Ganzi  and  Qiangtang  terranes,  we 
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observe  high  isotropic  1/QLg  (1/400-1/150  at  1  Hz).  A  low  1/QLg  anomaly  is  also  observed  near  the  Eastern 
Himalayan  Syntaxis,  which  is  probably  consistent  the  suggestion  of  an  underthrusting  cold  Indian  crust. 

CONCLUSIONS  AND  RECOMMENDATIONS 


In  this  project,  we  are  processing  a  large  set  of  new  broadband  waveform  data  to  obtain  three  velocity 
models  using  Pn,  body  waves,  and  surface  waves  as  well  as  a  crustal  Q  model  using  a  new  technique  for 
major  portions  of  western  China.  The  three  separate  velocity  models  show  consistent  results  including  a 
low  velocity  in  the  northern  Qiangtang  and  the  Kunlun  belt  and  a  high  velocity  anomaly  within  the  Qaidam 
Basin.  This  suggests  a  strong  lateral  variation  and  vertical  coherency  of  velocity  structures  for  the  whole 
lithosphere  in  eastern  Tibet  and  adjacent  regions.  A  high  seismic  attenuation  zone  is  observed  along  the 
Kunlun  belt  where  the  seismic  attenuation  has  a  strong  azimuthal  anisotropy.  The  isotropic  seismic  Q, 
suggested  as  an  estimate  of  crustal  Q,  is  low  in  northern  Qiangtang  and  the  Kunlun  belt  and  high  in  the 
Qaidam  Basin,  which  corresponds  with  the  crustal  velocity  models.  The  consistent  patterns  of  direction 
distribution  of  Rayleigh  wave  azimuthal  anisotropy  at  different  depth  and  Q  azimuthal  anisotropy  suggests 
a  coherency  of  deformation  between  the  crust  and  upper  mantle  in  this  region. 
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